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ABSTRACT. Bifurcation theory is used to analyze the nonlinear dynamic 
stability characteristics of an aircraft subject to single-degree-of- 
freedom pitching-motion perturbations about a large mean angle of attack. 
The requisite aerodynamic information in the equations of motion can be 
represented in a form equivalent to the response to finite-amplitude 
pitching oscillations about the mean angle of attack. It is shown how 
this information can be deduced from the case of infinitesimal-amplitude 
oscillations. The bifurcation theory analysis reveals that when the 
mean angle of attack is increased beyond a critical value at which the 
aerodynamic damping vanishes, new solutions representing finite- 
amplitude periodic motions bifurcate from the previously stable steady 
motion. The sign of a simple criterion, cast in terms of aerodynamic 
properties, determines whether the bifurcating solutions are stable 
(supercritical) or unstable (subcritical) . For flat-plate airfoils 
flying at supersonic/hypersonic speed, the bifurcation is subcritical. 



implying either that exchanges of stability between steady and periodic 
motion are accompanied by hysteresis phenomena, or that potentially large 
aperiodic departures from steady motion may develop. 



1. INTRODUCTION. Problems of aerodynamic stability of aircraft flying 
at small angles of attack have been studied extensively. With increasing 
angles of attack the problems become more complicated and typically 
involve nonlinear phenomena such as coupling between modes, amplitude 
and frequency effects, and hysteresis. The need for investigating sta- 
bility characteristics at high angles of attack was clearly demonstrated 
by Orlik-Ruckemann [1] in his survey paper which largely deals with 
experiments . 

On the theoretical side, the greater part of an extensive body of 
work is based on the linearized theory, in which the unsteady flow is 
regarded as a small perturbation of some known steady flow (possibly 
nonlinear in, e.g., the angle of attack) that prevails under certain 
flight conditions. The question of the validity and limitations of such 
a linearized perturbation theory is of fundamental importance and yet 
has been investigated only rarely. One may argue that in principle, it 
should be possible to advance to higher and higher angles of attack a 
by a series of linear perturbations , since the solution at each step 
should include a steady-state part which, when added to the previous 
steady-state solution, would provide the starting point for the next 
perturbation. This may well be true provided that at each step the 
steady motion is stable both statically and dynamically, and that the 
actual disturbances, e.g., the amplitude of oscillation, remain small. 
However, the linear perturbation procedure must eventually cease to be 
valid when the angle of attack exceeds a certain critical value o ^ 
at which the steady motion is no longer stable. 



In this paper we investigate the stability characteristics of an 

aircraft trimmed to a mean angle of attack a near a at which the 

m cr 

steady motion becomes unstable. Padfield [2] studied a similar problem, 
using the method of multiple scales, which is valid only for weakly non- 
linear oscillations. We shall study the problem by means of bifurcation 
theory. This will allow us to draw on recent mathematical developments 
(e.g., [3]) that are particularly well suited to investigating funda- 
mental questions in linear and nonlinear stability theory. A numerical 
scheme based on bifurcation theory was proposed earlier [4] for analyz- 
ing aircraft dynamic stability in a rather general framework. More 
recent work [5] demonstrates the considerable potential of bifurcation 
theory in flight dynamics studies, particularly toward establishing a 
method for the design of flight control systems to ensure protection 
against loss of control. On the other hand, while acknowledging the 
importance of the aerodynamic model in determining the aircraft stability 
characteristics, neither of these works contains an adequate assessment 
of the model requirements. The treatment of unsteady flow effects, in 
particular, receives no attention. In contrast, we shall focus on just 
this aspect of the problem at the expense of narrowing the scope of the 
motion analysis. Restricting the motion to a single-degree-of-freedom 
pitching oscillation will enable us to analyze a motion for which com- 
plete aerodynamic information is available, for certain aerodynamic 
shapes, in the form of exact solutions of the inviscid supersonic/ 
hypersonic unsteady flow theory [6-10]. In this way it will be possible 
to establish a form revealing a precise analytical relationship between 



the basic aerodynamic coefficients and the characteristics of the motion. 
For a more extensive account of the analysis presented here, see [11]. 

2. MATHEMATICAL FORMULATION. 

A. The Coupled Dynamic /Aerodynamic System. We restrict attention 
to the single-degree-of-freedom pitching oscillations of an aircraft 
about a fixed trim angle. Before time zero, let the aircraft be in 
level, steady flight with its longitudinal axis inclined from the hori- 
zontal velocity vector by the fixed trim angle of attack 0 m - At 
time zero the aircraft is perturbed from its trim position, but during 
the subsequent motion the center of gravity continues to follow a recti- 
linear path at constant velocity V . The instantaneous angle of attack 
a(t) is again measured relative to the horizontal velocity vector, while 
the inclination of a(t) from the fixed trim angle a , a(t) - a , will 
be designated £(t). The equations of motion are 

4- (o - a ) - 1 (2.1a) 

dt m 

I = M(t) (2.1b) 

where I is the moment of inertia and M(t) the instantaneous aerody- 
namic pitching moment, both referred to the center of gravity. We 
assume that the moment required to trim the aircraft at o m has been 
accounted for, so that M(t) is a measure of the perturbation moment 
only. 

Now we consider the equations governing the flow around the aircraft, 
assuming for simplicity that it is permissible to neglect viscous effects. 
The inviscid unsteady flow-field equations are 



(2.2a) 
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where p, p, and v are the pressure, the density, and the velocity 
vector, respectively, and y is the ratio of specific heats of the 
medium. Denote the position vector originating from the center of 
gravity by r and let the equation of the body surface be 
B(r,£(t)) = 0, which, as noted, depends on the instantaneous displace- 
ment angle 5(t). The tangency condition at the body surface is then 

+ v • VB = 0 on B - 0 (2.3) 


which also depends on £(t). The far-field boundary conditions (relative 
to a coordinate system fixed in the aircraft) are 

P = Poo’ P = Poo’ V - (2.4) 

for subsonic flight, or the well-known shock jump conditions for super- 
sonic flight. 

It is clear from the above that the pitching motion of the aircraft 
£(t) is coupled to the unsteady flow field (e.g., the pressure p) 
through (2.3), and, more directly, through M(t) in (2.1a&b), which is 
determined from the instantaneous surface pressure through 

S(t) - JJ ? * P dS (2.5) 

S:B=o 

The principal difficulty, however, is that the instantaneous surface 
pressure in (2.5) depends not only on the current state of . the flow 





field but also on its prior states. This means that past values of E, 
figure in the determination of the current state of the flow field and 
M(t) is a functional , not a function, of £(t). Thus, unless approxi- 
mations are introduced, a determination of E, starting from given 
initial conditions requires the simultaneous solution of the coupled 
equations (2. 1-2. 5). 

Although the simultaneous solution of the coupled equations (2. 1-2. 5) 
in principle represents an exact approach to the problem of determining 
time-histories of maneuvers from given initial conditions, it is 
inevitably a difficult and costly approach (see the discussion of [12]). 
Approximate approaches leading to simpler, less costly computations are 
a practical necessity. To date, these computations have invoked the 
assumption of a slowly varying motion of known form (e.g., a harmonic 
motion) whose aerodynamic force and moment response at t could be 
calculated. This stratagem, in effect, uncouples the flow-field equa- 
tions from the equations of motion. In a series of papers (see in 
particular [13]), Tobak and Schiff have shown how this stratagem can be 
rationalized to create mathematical models of the aerodynamic response 
at various levels of approximation of the dependence on the past motion. 
(See also the contribution of Tobak, Chapman, and Schiff in this 
collection.) 

B. Uncoupling Near Neutral Dynamic Stability Boundary. Let us 
suppose that there is a mean angle of attack = a ^ at which the 
damping- in-pitch coefficient vanishes. Then, at and beyond this angle 
of attack, the steady motion of the aircraft will no longer be stable 


so that, in response to a disturbance, the motion will seek a new stable 
state which usually will consist of a finite-amplitude periodic oscilla- 
tion about the mean angle of attack. The changeover from a stable 
steady motion to a stable periodic motion is called a Hopf bifurcation 
(see [3, 14], also Sec. 4). Thus, in the vicinity of a cr » where the 
damping- in-pitch coefficient will be near zero, the consequent persis- 
tence of oscillatory motions will make it particularly appropriate to 
assume a periodic past motion about a mean position in calculating the 
aerodynamic pitching moment from the flow-field equations for use in 
(2.1). This assumption is consistent with the mathematical modeling 
approach of Tobak and Schiff [13] at the second level of approximation. 

It will be convenient to write M(t) in the form 


M(t) - -y p V 2 S£[C (§,i,0 ) 

2 °° oo m m 


- C (0,0, a )] 
m ’ * m /J 


( 2 . 6 ) 


The function 0^(5, £,a m ) is the pitching-moment coefficient resulting 
from a finite-amplitude periodic pitching motion £(t) about the center 

of gravity, and 0^(0, 0,a m ) is its steady-state value at the mean angle 

• 

of attack a . The function C (£,£,a ) depends on the instantaneous 
m m m r 

• ' 

displacement angle ?(t), its rate of change £(t), and the mean angle 

of attack a . It also depends, of course, on the location h of the 
m 

center of gravity relative, say, to the wing leading edge, the flight 
Mach number M , the ratio of the specific heats of the medium y, and 
the aircraft shape. A great amount of work has been devoted to the 
theoretical determination of the function for the case of pitching 

oscillations of infinitesimal amplitude. For example, for the cases of 
a wedge, a flat plate in supersonic /hypersonic flow [7, 8] and an airfoil 
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of arbitrary profile in the Newtonian limit [9], this function is 

available in exact analytical form for large, as well as small, mean 

angles of attack a , and includes the critical angle a . On the 

m cr 

other hand, little is known about the pitching-moment coefficient C 

m 

when the amplitude of oscillation is finite, except for the special 
case of a slowly oscillating wedge [10] . In the next section we shall 
show how the function for slow pitching oscillations of finite 

amplitude may be obtained from its behavior in the limit of oscillations 
of infinitesimal amplitude. 


C. The Pitching-Moment Function for Slow Pitching Oscillations 
of Finite Amplitude. Consider a uniform flow V past an aircraft that 
is undergoing a slow pitching oscillation of finite amplitude about a 
mean angle of attack a . The instantaneous angular displacement from 
the mean position is measured by £(t) so that [ g [ mav is the finite 
amplitude of the oscillation. It is required to calculate the unsteady 
pressure field, hence the form of the pitching moment M(t) in (2.1). 

Let the cylindrical coordinates (r,<|>,z) be such that the z-axis coin- 
cides with the lateral axis through the center of gravity. Let the 
equation of the body surface at the mean angle of attack cr^ be 
<f) = Cm + A(r,z). Then its equation at time t is 

B(r,<j>,z,t) = A(r,z) + a m + £(t) - <j> = 0 (2.7) 

With velocity vector v = ( v r > v ( j ) > v z ) > the boundary condition (2.3) 
becomes 


v , 


- r (5 


+ v 


3A 

9r 


+ v 


( 2 . 8 ) 
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at <j> = A(r,z) + a m + 5(t). Equations (2.2), (2,4), and (2.8) complete 
the formulation for calculating the flow field in terms of 5(t). Now, 
for a long-established periodic motion of the aircraft, the aerodynamic 
response is also periodic. Furthermore, invoking the assumption of 
slow oscillations allows us to neglect terms of 0(5 ,5) and higher and 
to write 

• • . 

P = P 0 + £Pi > p = p 0 + 5p x , v = v Q + 5v x (2.9) 

where ( ) 0 and ( ) x are independent of |(t), but may depend on time 
through the function 5(t). Since the zeroth-order quantities ( ) 0 
represent the flow field when 5 = 0, such a flow must be quasi-steady; 
i.e., the time variable t can only appear implicitly through the 
instantaneous displacement angle 5(t). Hence, 

P 0 = P 0 (a m + 5(t) ,r,<(),z) , etc. (2.10) 

To derive the mathematical problem for the first-order quantities 

( ) l5 we substitute (2.9) into (2.2), (2.4), and (2.8), use (2.10), and 

• « •• 

neglect terms of 0(5 ,5). Thus, 


V • (p„v + p vj 
i K i o 
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(2.11a) 

(2.11b) 


(2.11c) 
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and 



•bz) 


(2. lid) 


(2. lie) 


at 

<j> = A(r,z) + a + ?(t) 
m 

Since the time variable t appears in (2.11) only as a parameter 
through 5(t), the solution to (2.11) must be of the form 

p i = p i ( a m + » etc - (2.12) 

The forms of p Q and p^ in (2.10) and (2.12) determine the form of the 
pressure p in (2.9) which, in turn, determines the form of the 
pitching-moment coefficient after the integrations indicated in (2.5) 
have been carried out. Thus 

C a, La) - f (cr + O + & g(a + O (2.13a) 

m m m V m 

00 

In the case of oscillations of infinitesimal amplitude, (2.13a) reduces 
to 

C - f(o ) + Zf'(o ) + g(o ) (2.13b) 

m m m V m 

CO 

It is now clear that the functions f’(a ) and g(a ) are related to the 

stiffness derivative ^(a^) and the damping-in-pitch derivative D(o m ) 

at an angle of attack a , as defined in classical aerodynamics, by 

= -s ( a m ) » S(° m ) = -D(a) (2.14) 

m m m m 

Comparing (2.13a) and (2.13b), we conclude that knowing the stiffness and 
damping derivatives from the results of calculations for oscillations of 
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infinitesimal amplitude enables one to obtain immediately the pitching- 
moment coefficient for the corresponding finite-amplitude case. 

This general conclusion is supported by the form of the exact solution 
for C m for the wedge oscillating at large amplitude given in [10], 

(To correct a misprint in [10], a term hA. [cos(0 - 0 ) - 1] should be 

** 3 

added to the right-hand side of Eq. (12b) in [10].) 

Having determined an appropriate form of M(t) for use in the iner- 
tial equations of motion (2.1), we now rewrite the equations introducing 


(2.13a) along with dimensionless time (i.e., characteristic lengths of 
travel) t = V^t/A. (Note that we shall retain the symbol ( ) to desig- 
nate a time derivative. Henceforth, however, the derivative will be 
with respect to dimensionless time: ( ) = d/dx.) Let 


F(5,5,0 = - ■ v - ■ ' ■■ ■ = k[C (5,5,0 - C (0,0,o ) ] 
m y £) 2 m m m m 


= k [ f (a + 5) - f (a ) + 5g (a +5)] (2.15) 

m m m 


with 


p Sit 3 
00 

21 * 



The inertial equations of motion become 

ff = 1 (2.16a) 

ff- = F(5,5,a m ) (2.16b) 

• • 

An expansion of F(5,5,O m ) in a Taylor series in 5 and 5 and a 

change of notation ^ “ 5-» u “ 5 yields for (2.16) 
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u, = A. . (a )u. + B. (a )u.u. + C. M (a )u.u, u. 
i ij m 2 ljk m j k ijk£ m j k £ 


+ 0 ( ul 4 ) , (1 = 1,2) 


where 


A = 


,-KS(cr ) -kD(o ), 
m m 


B ijk ~ 0 5 


B„ = 


1 3 Z F 


zjk 2! 3^91^ 


u=o 


(2.17) 


(2.18a) 


(2.18b) 


C ljk£ ° * C 2jk£ 


1 

3! 


3 3 F 

9u j 3u k 9u £ 


->• 

u=o 


(2.18c) 


(Although (2.16a&b) have been derived on the assumption of slow oscilla- 
tions (terms in of 0 (| 2 ,?) omitted), our subsequent bifurcation 

analysis of (2.16) will hold for general F(£,£,c m ), i.e., as if no 
restriction had been placed on the magnitude of £.) 

In (2.18) the tensors B and C represent the effects of finite ampli- 
tude to the second and third order. We note that the following symmetry 
properties hold: 


B 2jk B 2kj 

C 2jk£ = C 2j£k " C 2kj£ = C 2k£j = C 2£jk = C 2£kj 
On the basis of (2.17), we shall study the linear and 
of the motion in subsequent sections. 


(2.19a) 

(2.19b) 

nonlinear stability 


3. LINEAR STABILITY THEORY. The stability of the steady motion at 

an angle of attack a to infinitesimal disturbances is determined by 

m 

the nature of the eigenvalues of the matrix A. They are 



(3.1) 


Ai(a ) = -y [-kD(o ) ± /k 2 D 2 (cr ) - 4«S(a )] 

2 m / m m m 

Case I : S(cr m ) <0. In this case A 1 >0, X 2 <0. The steady motion 

at this angle of attack 0 m is always unstable. 

Case II: S(o ) > 0. 

m 

Case Ila : D(o ) < 0. In this case Re (A.) > 0 and the steady motion 
m i J 

at a is unstable, 

m 

Case lib: D(a ) > 0. In this case Re(Xi) < 0 and the steady motion 
m 2 

at a is stable, 

m 

Thus, only in Case lib, when both stiffness and damping-in-pitch 

derivatives are positive, is the steady motion at angle of attack 

stable to infinitesimal disturbances. In fact, stability theory [3] can 

be used to show that stability of the steady motion in this case is 

assured only if the disturbance is sufficiently small. 

In all cases in which the linearized theory predicts growth of the 

disturbance amplitude, the growth predicted is of exponential form and 

hence must cease to be valid after some finite time when the amplitude 

is no longer small. Thus, what eventually happens to a motion for which 

linearized stability theory predicts an initial growth of disturbances 

cannot be determined from the linearized theory itself. Instead, the 

full nonlinear inertial equations of motion, or a suitable approximation 

of them, such as (2.16), must be adopted to determine the ultimate state 

of the motion. Of particular interest is the dynamic stability boundary 

a = a , where S(a ) > 0, D(cr )=0. The stability characteristics 
m cr cr cr' J 

near this boundary will be studied in the next section. 


4. NONLINEAR STABILITY THEORY. At the dynamic stability boundary 

a = a , we have S(a ) > 0 and D(a ) = 0, hence 
cr cr' cr 

Xi(a ) = ±i/icS(a ) = ±io) 0 (4.1) 

2 Cl. Cl w 

The existence of purely imaginary eigenvalues of the matrix A at 

o = o is the characteristic sign of a Hopf bifurcation [3, 14], 

signaling a changeover from stable steady motion to periodic motion. 

On crossing a m = cr , t * ie stea ^Y m °tion that had been stable for 

a < a will become unstable to disturbances, resulting (after a 
m cr 

transient motion has died away) in the existence of a new motion, which 
(if it is stable) will be periodic. In the vicinity of a m = o , the 
circular frequency of the periodic motion will be nearly equal to oj 0 . 
We call the new solution of the equations of motion a bifurcation solu- 
tion. In this section we shall determine its character and a criterion 
for its stability. 

For slightly larger than a cr > the eigenvalues of the matrix A 

are 

Xi = - \ kD (0 ± (4.2) 

2 2 m m 

where 

fi(a ) = /kS(o ) - k 2 D 2 (cj )/4 (4.3) 

m m m 


We shall assume that 



which is the usual case in applications [6]. (The case D’(cr ) > ® can 
be treated in exactly the same way.) The normalized eigenvector "£(a ) 
associated with the eigenvalue A(o m ) is 


t(a) = 

m 


<or_> 


m 


1 - i 


Z 2 (a) l/uiTT U (a) 


(4.5) 


m 


m' 


m 


whereas the adjoint eigenvector with eigenvalue ^(cr^), which is 


the complex conjugate of 1 ( 0 ^), i s 


t*(o) 

m 


V^V 


1 + i 

2/Q(a ) 

m 


X(a ) + kD(ct ) 
m m 


(4.6) 


A. Hopf Bifurcation. The bifurcation solution u ( T > a m ) may be 


written as 


u = a(r)t + a(t)t 

Following Iooss and Joseph ([3], p. 125), we get 
a = (s) + e 2 b 2 (s) + e 3 b 3 (s) 4- 0(e 1 *) 

s = [w. + e 2 w, + 0(e 1 *)]x 


a = a + e 2 o + 0(e H ) 
m cr m 2 v 


(4.7) 


(4.8) 


where, for brevity, we omit the lengthy solution forms for e, b^ , cu 2 . 


n 


and o m2 (cf. [11]). The solution is periodic in t with circular 


frequency equal to u) Q + e 2 w 2 + 0(e 4 ), 


B. Stability of the Periodic Bifurcation Solution. According to 
Floquet theory [3], the stability of the periodic bifurcation solution 
(4.8) is determined by the sign of an index y. To 0(e It ), y has the 


form 


and the periodic bifurcation solution is stable if y < 0, unstable if 


y > 0. Since we have assumed D'Ca^) < 0, stability thus depends on 
the sign of o m2 , with <J m2 > 0 denoting stability and o m2 < 0 insta- 
bility. It remains to cast o m2 in more recognizable terms. After 
considerable manipulation, we get 

i 3i 


y = — 


7l!l 4. 2 3 2 F \ 3 2 F 2 / 3 3 F 

\3u 2 0 3u 2 / 3u x 3u 2 0 ^3u 2 3u 2 


F 2 3 F 

+ or r 

0 3u; 


)]* 

/J u=o 


a -a 
m cr 


(4.10) 

in terms of F( u i ,u 2 ,a m ) . From (2.15) we see that the function F is 
directly related to the pitching-moment coefficient 0^(5, 5 , 0 ^) acting 
on the aircraft which is performing a finite-amplitude pitching oscilla- 
tion 5 around a mean angle of attack o . Equation (4.10) demonstrates 
that the stability of the periodic motion near the dynamic stability 
boundary a is determined by the behavior of the aerodynamic response 

C (£,5,a ) in that vicinity, 
m m 

With the assumption of slow oscillations under which the form of 

4 a tl 

(2.15) was derived (terms of 0(5 ,0 neglected), we may substitute 
(2.15) into (4.10) to get 

_ *2,A 

(4.11) 


e 2 / 2 3 2 f 3g 2 3 V 

4(og V 35 0 H / 5=0 




a =o 

m cr 


From (2.14) and the structure of the functions f and g, we write 
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9f 

95 



= -S (a ) , 

9 2 f 


= -S ' (a ) 

5“u 

m 

95 2 

5=0 

m 

5=o 

- -D(a m ) , 

9g 

95 

5=o 

= -D’ (a ) , 
m 


1 


lis 

95 2 


- -D"(cr ) 
5=0 m J 


(4.12) 


By using (4.2) and (4.12) in (4.11) (with S(a cr ) > 0), we arrive at the 


simple form 


y 


-e 2 (kS ) 1 / 2 
4 




(4.13) 


m cr 

(The simplicity of this result suggests that it might be possible to 
derive it by a less formal method. This is indeed the case: We have 

verified the result by a physical approach familiar to workers in the 
field of nonlinear mechanics.) Equation (4.13) reveals that the sign of 
y, and thus the stability of the bifurcation solution, is independent of 
the scalar and inertial properties of the aircraft. Rather, stability 
depends only on whether the aerodynamic property (D'/S) is increasing 
or decreasing on crossing the dynamic stability boundary = .a . ■ The 
two possibilities are well illustrated in the form of bifurcation dia- 
grams as shown in Figs, la and lb. 

In a bifurcation diagram, the abscissa represents the parameter that 
is being varied, in our case the mean angle of attack a . The ordinate 
is a parameter characteristic of the bifurcation solution alone. In our 
case it is e, a measure of the amplitude of the periodic bifurcation 
solution. Stable solutions are indicated by solid lines, unstable solu- 
tions by dashed lines. Thus, over the range of mean angle of attack 

a <a where the steady-state motion is stable, e is zero, and the 
m cr J 



stable steady motion is represented along the abscissa by a solid line. 

The steady motion becomes unstable for all values of a > a as the 

m cr 

dashed line along the abscissa indicates. Periodic solutions bifurcate 

from a m = a ^ either supercritically or subcritically. 

When (d/da ) (D '/S) | _ >0 (implying a m2 > 0), the bifurcation 

is called supercritical and its characteristic form is shown in Fig. la. 

Stable periodic solutions (solid curves in Fig. la) exist for values of 

a m - a cr > 0. The amplitude of the periodic solution at a given value 

of a m - a is proportional to e, hence is vanishingly small when 

°m ~ °cr sma ^> varying essentially as (a ffl - a^) 1 / 2 . 

When (d/d a )(D'/S)| _ <0 (implying a mz < 0), the bifurcation 

m a m -a cr 

is called subcritical and its characteristic form is shown in Fig. lb. 

Periodic solutions exist for values of a -a <0, but they are 

m cr 

unstable (dashed curve in Fig. lb). Whether stable periodic solutions 

do or do not exist for a m > a cr depends predominantly on the behavior 

of the damping- in-pitch derivative D(a ) for a > a . If no such 
r r m m cr 

stable periodic solutions exist for a > a , then when the mean angle 

m cr 

of attack a m is increased beyond a cr the aircraft may undergo an 
aperiodic motion whose departure from the steady motion at a^ = a cr 
is potentially large. 

In the more likely event that stable periodic solutions do exist for 

CT m > °cr ^ an exam Pl e gi ven later), their amplitudes must be finite, 

and not infinitesimally small, even for small positive values of 

a -a . It is likely that this branch of stable periodic solutions 
m cr 

will join that of the unstable branch in the way illustrated in Fig. lb. 



In this event, the form of the bifurcation curve for values of a < a 

m cr 

helps explain the situation mentioned earlier, where it was noted that 
the steady-state motion could be stable to sufficiently small distur- 
bances but become unstable to larger disturbances. Thus, Fig. lb 
suggests that for cr^ < cr^, so long as disturbances are of small 
enough amplitude to lie below those of the unstable branch of periodic 
solutions (curve OB in Fig. lb), they will die out with time and the 
steady motion will remain stable. However, disturbances with ampli- 
tudes sufficiently larger than those of the unstable branch may actually 
grow up to the ultimate motion as x, •* «>, which will be that of the 
stable branch of periodic solutions (curve BA in Fig. lb). Finally, we 
note that if the motion does attain the stable branch of periodic solu- 
tions (say, for a < a ) then hysteresis effects will manifest them- 
m cr 

selves with further changes in a . When a is increased beyond a , 

e m m 1 cr 

the motion will continue to be periodic with finite amplitude (point A 
in Fig. lb). If o m is now decreased below cr , the periodic motion 
will persist, even at values of a m where previously there had been 
steady motion when a was being increased. Not until a is dimin- 
ished beyond a certain point (point B in Fig. lb) will the motion return 
to the steady-state condition (point C in Fig. lb) that had been experi- 
enced when a m was increasing. 

To further explore the implications of (4.13), we invoke some approx- 
imate relationships between the damping- in-pitch derivative D(cr m ) an( l 

the stiffness derivative S(o ). Tobak and Schiff have argued (see, in 

m ■ 

particular, [15]) that, to good accuracy, a linear relationship should 
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exist between D and S at any value of a ; i.e., D(a ) * a - bS(a ), 

J m’ m m y 

with b > 0. In addition, we require that D vanish at a = a 

m cr 


The two requirements yield the form 

D(a ) = b[S(a ) - S(a )] , b > 0 
m' v cr y m 


(4.14) 


(The validity of the approximate relation (4.14) has been verified in 
the present study (via comparison with exact results in [7]) for use 
with oscillating flat-plate airfoils in supersonic /hypersonic flow and 
in [12] for use with oscillating flaps in transonic flow.) Replacing 
D in (4.13) by (4.14) casts the criterion solely in terms of S: 

_ e 2 (<S) 1 / 2 b d 2 


V = 


da" 


[In S(a ) ] 


m 


(4.15) 


m 


a -o 
m cr 


Thus, near a = a , if the stiffness derivative S(a ) increases with 
m cr m 

a^ slower/faster than exponential, the periodic bifurcation solution 
is stable (supercritical) /unstable (subcritical) . Cast in terms of D 
instead of S, the criterion states that decreasing D with respect to 
slower /faster than exponential results in stable (supercritical)/ 
unstable (subcritical) periodic bifurcation solutions. 

Examples of ^(a^) and D((J m ) that obey (4.14) and exhibit the various 
possibilities are as follows. 

Example 1 : 

S i (a m ) = S o ex P£ k < a m “ a 

D (o ) = bS {1 - exp[k(a 
i m O' K m cr 


) + m(a - a ) 2 ] ”) 

m cr I 

-a ) + m(a - a ) 2 ]}J 

cr m cr J 


(4.16) 


Example 2 : 

S 0 (a ) = S. exp[k(a - a ) + m(a - a ) 2 - n(a — a )**] ^ 

2 m 0 ' m cr m cr m cr l 

D,(cr ) = bs o^ “ e x Ptk(cr “ a ) + "»(a - a ) 2 - n(a - a ) 1 *]} J 

<4 m u m cr m cr m cr J 

(4.17) 

where S Q , k, and n are positive constants. According to (4.15), 
m > 0 corresponds to unstable periodic bifurcation solutions (sub- 
critical bifurcation), and m < 0 corresponds to stable periodic bifur- 
cation solutions (supercritical bifurcation) . In the case of subcritical 
bifurcation (m > 0) in Example 1, the damping-in-pitch derivative Dj (cr ) 
continues to decrease to very large negative values as (o^ - ° cr ) 
increases, which makes it very unlikely that the system (2.16) will 

have a stable periodic motion as a solution for a > a . On the other 

m cr 

hand, in Example 2 the damping-in-pitch derivative D. (a ) becomes posi- 

^ m 

tive for sufficiently large | o m - a cr | so that a stable periodic 
motion may be a possible solution for a > a , m > 0. Indeed, S. 
and D 2 in Example 2 fulfill all of the conditions set by a theorem of 
Filippov [16], the satisfaction of which guarantees the existence of a 
stable periodic motion of the system (2.16) for > ° cr « As we have 
noted, resulting as it does from a subcritical bifurcation, the solution 
will exhibit hysteresis effects with variations in a below a 

5. SUPERSONIC /HYPERSONIC FLAT-PLATE AIRFOILS. To illustrate the 
application of bifurcation theory in a concrete case, we consider a 
flat-plate airfoil in supersonic /hypersonic flow. The stiffness deriva- 
tive S(cr ) and the damping-in-pitch derivative D(o ) are known as 
m m 
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analytic functions of o^ up to the angle-for-shock detachment [6-8]. 

They depend on the flight Mach number M^, the ratio of specific heats 

y (here taken to be that of air, y = 1.4), and the (dimensionless) 

distance h of the center of gravity from the leading edge, here taken 

as a fraction of the chord length &. Results are presented in Fig. 2 

of log S(o ) versus a (with M = 2.0, h = 0) and in Table 1 of the 

index y for various combinations of and h. It is shown in Fig. 2 

that S(o ) increases faster than exponential near a = a = 15.75°, 
m r m cr 

and in Table 1 that y is always positive. We conclude that whenever 
the flat-plate airfoil becomes dynamically unstable [D(c? cr ) = 0], the 
ensuing bifurcation always will be subcritical. 

Accordingly, there are two possibilities. One, a subcritical bifur- 
cation curve such as that sketched in Fig. lb exists, in which case the 

airfoil motion will find stability at values of a > a in a 

m cr 

periodic oscillation of finite amplitude. At values of a m < a , the 

steady-state motion will be stable for small disturbances, but for 

larger disturbances the airfoil motion will again seek the stable 

periodic motion. Exchanges between these two stable modes at a < a 

m cr 

will be accompanied by hysteresis effects. Two, alternatively, a bifur- 
cation curve such as that sketched in Fig. lb does not exist, in which 
case a potentially large aperiodic departure from the steady-state 
motion may occur as exceeds cr , In either case, on exceeding 

a cr the loss of stability of the steady-state motion must entail a 
discrete change to a new stable condition. 


# 


6. CONCLUDING REMARKS. Bifurcation theory has been used to analyze 
the nonlinear dynamic stability characteristics of an aircraft subject 
to single-degree-of-freedom pitching-motion perturbations about a large 
mean angle of attack. Setting up the equations to which the bifurcation 
theory was applied required, first, determining conditions under which the 
inertial equations of motion and the gas-dynamic equations governing the 
flow could be decoupled and, second, showing how the required aerody- 
namic responses to finite-amplitude oscillations could be obtained from 
the responses to infinitesimal-amplitude oscillations. 

Results of the bifurcation theory analysis revealed that when the 
mean angle of attack is increased past the critical point where the 
aerodynamic damping vanishes, new solutions describing finite-amplitude 
periodic motions bifurcate from the previously stable steady motion. 

The sign of a simple criterion, cast in terms of aerodynamic properties, 
determines whether the bifurcating solutions are stable (supercritical) 
or unstable (subcritical) . For flat-plate airfoils flying at supersonic/ 
hypersonic speed, the bifurcation is subcritical, implying either that 
exchanges of stability between steady and periodic motions will be 
accompanied by hysteresis phenomena, or that a potentially dangerous 
aperiodic motion may develop. In either case the loss of stability of 
the steady-state motion must be accompanied by a discrete change to a 


new stable state. 
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Figure Captions 

Fig. 1 Typical forms of bifurcation diagrams near the dynamic stability 
boundary O cr where D ( a cr ) = 0. 

(a) Supercritical, (d/dcr^) [D* ( G m >/S( 0 m ) l a =a > 0. 

(b) Subcritical, (d/da )[D'(a )/S(o )] < 0. 

m m m o m -<T cr 

Fig. 2 Stiffness derivative S(a ) versus mean angle of attack a for 
° m m 

a flat-plate airfoil in a supersonic free stream: = 2.0, h = 0, 

Y = 1.4, K = 1. 
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Table 1 Values of stability criterion h) for flat-plate airfoil: 

k = e 2 =1, y = 1.4; u > 0 subcritical bifurcation, y < 0 supercritical 
bifurcation 
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